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Abstract 

Using the quark-meson coupling (QMC) model we study nuclear matter 
from the point of view of quark degrees of freedom. As the nucleon model 
we adopt the MIT bag model and the relativistic constituent quark model, 
where a square well and harmonic oscillator potentials are used to confine 
the quarks. We introduce the Lorentz-vector type confining potential as well 
as the Lorentz-scalar type one in order to examine how the vector confining 
potential contributes to the properties of the nucleon and nuclear matter. 
Next, we perform a re-definition of the scalar field in matter and transform 
the QMC model to a QHD-type model with a non-linear scalar potential. The 
result obtained from QMC is then compared with the potentials which are 
determined so as to fit various properties of finite nuclei and nuclear matter 
in relativistic mean-field models. The QMC model provides the parameters 
k ~ 20 — 40 (fm _1 ) and A ~ 80 — 400 for the standard, non-linear scalar 
potential. We discuss a relationship between the QMC and QHD-type models 
in detail. 

PACS numbers: 24.85. +p, 12.39.Ki, 24.10.Jv, 21.65. +f 
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I. INTRODUCTION 



In the conventional nuclear physics we have assumed that the nucleon properties in the 
nuclear medium are not changed from those of the free nucleon. Thus, low- and medium- 
energy nuclear physics phenomenology has been successfully described in terms of (point- 
like) nucleons. However, after the EMC effect in the nuclear structure functions was ob- 
served [0 there has been considerable interest in the modification of hadron properties in 
the nuclear medium. It is one of the central questions in nuclear physics how the effect of 
the quark substructure of hadrons emerges in nuclear phenomena. 

At present rigorous studies of quantum chromodynamics (QCD) are limited to mat- 
ter system with high temperature and zero baryon density. Because of the complex non- 
perturbative features of QCD it is very difficult to calculate properties of nuclei at low- 
and medium-energy. Thus, we need to build models which help to bridge the discrepancy 
between nuclear phenomenology and the underlying theory. Such models may be necessarily 
crude because the nuclear many-body problem at the fundamental QCD level is intractable. 
However, it is very important and intersting to challenge this problem. Then, much effort 
has been devoted to the study of effective models for nuclear matter and finite nuclei based 
on various quark models 0. 

About a decade ago, Guichon || proposed a relativistic quark model for nuclear matter, 
where it consists of non-overlapping nucleon bags bound by the self-consistent exchange 
of isoscalar, scalar (a) and vector (u) mesons in mean- field approximation (MFA). This 
model has been further developed as the quark- meson coupling (QMC) model It was 

shown that the model can reproduce the saturation properties of nuclear matter and that 
it provides a fair description of finite nuclei (for recent reviews, see Ref. ||). Since in QMC 
the mesons couple not to the nucleons but to the quarks directly the nucleon properties are 
modified self-consistently, depending on the scalar field in a medium. 

On the other hand, recent theoretical studies show that various properties of finite nuclei 
can be described by Quantum Hadrodynamics (QHD) [|7| very well. Although the original 
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version of QHD [|7| provides the large reduction of the nucleon mass in nuclear matter, which 
is the source of the large spin-orbit splitting in finite nuclei, the model also gives the large 
nuclear compressibility around 550 MeV. In order to improve this Boguta and Bodmer || 
first introduced self-interaction terms for the a meson. It is now well established that such 
non-linear interactions are necessary to reproduce the bulk properties of finite nulcei and 
nuclear matter [f|. 

Our purpose here is to study a relation between QHD-type mean-field models and the 
QMC model and to examine how the internal structure of the nucleon sheds its effect on 
effective nuclear models. In QMC the nucleon mass in matter is given by a decreasing 
function of the a field. By applying a field re-definition to the scalar field in matter, the 
QMC model can be cast in a form equivalent to a QHD-type model with a non-linear scalar 



potential and a non-linear coupling to the gradient of the scalar field |T(J. Therefore, in 
principle, if the effective nucleon mass is known as a function of the scalar field, it could 
provide a prediction of non-linear interactions of the scalar field in nuclear matter through 
the procedure of the field re- definition. Inversely, it may imply that if a form of the scalar 
self-interactions in matter could be determined by analysing various experimental data we 
could get significant information on the quark substructure of the in-medium nucleon. 

The QMC model relies on the choice of a quark model of the nucleon. The MIT bag 
model has been used as the nucleon model in our previous works [|3|-f|. In the bag model 
the inside of the bag is assumed to be the perturbative vacuum and the quark moves with a 
current quark mass, which is nearly zero. There is an alternative idea for the nucleon model, 
namely the constituent quark model [fTI |. In this model, the quark takes a constituent quark 



mass of about 300 MeV in the nucleon, which is offered by spontaneous chiral symmetry 
breaking. Recently, Shen and Toki |I2] have proposed a new version of the QMC model, 
where the constituent quark model with a harmonic oscillator potential was used, and studied 
nuclear matter and finite nuclei — they refer the model as the quark mean-field (QMF) 
model. Although the quark acquires a heavy mass in this picture it is still important to 
treat it relativistically because the quark kinetic energy is not small compared with the 



constituent quark mass. Hence, in this paper we will use both the bag model and the 
relativistic constituent quark model in QMC and study the properties of nuclear matter and 
self-interactions of the scalar field. 

The outline of the present paper is as follows: In Sec. [nj we present relativistic quark 
models for the nucleon, which can be solved analytically. Using those quark models we 



calculate the nucleon properties in Sec. |T|. We also present our numerical results for nuclear 
matter. In Sec. [IV], by performing a re-definition of the scalar field we transform the QMC 
model to a QHD-type model. Then, we discuss non-linear scalar potentials generated from 
the QMC model. In Sec. [V| we give some discussions and conclusion. 



II. QUARK MODELS FOR THE NUCLEON 

In this section we summarize relativistic quark models to describe the nucleon. First 
we discuss the constituent quark model with a square well (SW) potential or a harmonic 
oscillator (HO) potential to confine the quarks. We introduce a Lorentz- vector type potential 
as well as a scalar type one to examine the role of the vector type confining potential in 
nuclear matter. In the constituent quark model it is assumed that the quark moves in the 
Nambu-Goldstone vacuum and hence the quark takes a constituent quark mass (around 300 
MeV). Next, we review the MIT bag model shortly. In the bag model, the quark moves in 
the perturbative vacuum (or the Wigner phase) and it takes a small, current quark mass. 



A. Constituent quark model 

We consider a light (u or d) quark (of mass m q ) moving under a potential V(r). The 
Dirac equation for the quark field ip q (r) is then given by 

[?7 • d-m q -V(r)]ip q (r) = 0, (1) 

where we take 

V(r) = (l + p l0 )U(r), (2) 



with a confining potential U(r) and a parameter (3{> 0) to control the strength of the 
Lorentz- vector type potential. We assume that the shape of the Lorentz-scalar type confining 
potential is the same as that of the vector type one. 

We first consider a square well potential and find a solution for the quark field a la 
Bogolioubov |T3 |. The potential is then given by 



f for r < R 

U{r) = \M, for r>R, (3) 

where R is the radius of the spherical well and M is the height of the potential outside the 
nucleon. After finishing all calculations, we take the limit M — > oo. 

The solution for the quark field can be found easily ]T3 |. Inside the well the lowest energy 
state is 

Ur) = JL e ->Ejm \ (4) 



'Arc \i<J ■ rg{r) 

with the normalization constant N, spin state Xs and the quark energy E. Here the upper 
f(r) and lower g(r) components of the quark wave function are respectively given by the 
spherical Bessel functions, jo(r) and ji(r). (For details, see Appendix.) Outside the nucleon 
volume we find a similar form of the solution with modified spherical Bessel functions. 



Note that since the quark wave function is proportional to e ~ M V 1 ~l 32r outside the well the 
condition (3 < 1 is necessary to confine the quark. 

We then demand that / and g must be continuous at r = R. This gives an eigenvalue 
problem, and we obtain a matching condition in the limit M — > oo: 



[l-{3)(E-m q ) 



3i(x), (5) 



\ (l + P)(E + m q y 

with the eigenvalue of the confined quark, x. The total energy is given by 3a/ R, where 
a 2 = x 2 + A 2 and A = Rm q . As in the MIT bag model there exist the spurious center of 
mass (cm.) motion and gluon fluctuation corrections etc. Here we add the familiar form, 
—z/R, with a parameter z for those corrections to the energy [13]. Thus, the nucleon mass 
(at rest) is given by 



Using the solution Eq.([|) we can calculate various nucleon properties analytically (see Ap- 
pendix) . 

This system may be described by a Lagrangian density 

C sw = Mil • d - m q )^ q d{R - r) - + p 1 ■ a)^ q 5{r - R), (7) 

where a M is the unit vector in time direction: a M = (1,0). This Lagrangian provides a 
boundary condition at r = R 

ij ■ ntpq = (1 + /fy • a)ip q , (8) 

where n M is the unit normal outward from the potential surface. This is equivalent to the 
condition Eq.([|) and ensures that the quark is confined permanently. (Note that actually 
the quark current flows in the azimuthal direction ||T4|| .) 

In the case where the confining potential is a harmonic oscillator potential with (3 = 1 
and c the oscillator strength 

V» = ^(l + 7o)cr 2 , (9) 

one can again find the quark wave function analytically. The solution to the Dirac equation 
is expressed by Eq.([|), where the upper and lower components are given in terms of gaussian 
functions [15|]. (For details, see Appendix.) The condition to determine the quark energy E 
is then obtained 



E + m q (E - m q ) = 3^c. (10) 

In the HO model the cm. energy can be estimated exactly as in the non-relativistic 
harmonic oscillator. It is just one third of the total energy |TjJ. Thus, the nucleon mass is 

M N = IE - E s , (11) 
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where E s describes gluon fluctuation corrections etc. Note that an equally weighted 
(Lorentz) scalar-vector potential does not produce a spin-orbit splitting in the baryon spec- 
trum, which may be consistent with the fact that the spin-orbit splitting in the observed 
baryon spectra is very small. 



B. MIT bag model 

In the MIT bag model the inside of the nucleon bag is considered to be the perturbative 
vacuum, and the quark takes a current quark mass. It requires extra energy to make a 
nucleon bag in which the quark moves freely, so we have to add a latent heat term to the 
Lagrangian density Eq.(0) 

C B = £sw-B9(R-r), (12) 

where the bag constant B describes the energy gap between the inside and outside of the 
nucleon. 

The quark wave function and the eigenvalue condition are given by the same forms as in 
the SW model. Thus, the nucleon mass is given as 

M N = ^ + £ B #, (13) 

where the z parameter is again introduced to take into account the sum of the cm. and 
gluon fluctuation corrections etc. 



III. QMC MODEL FOR NUCLEAR MATTER 

A. Effect of nucleon structure in a nuclear medium 

Here we consider how Eq.(|l]) is modified when the nucleon is bound in static, uniformly 
distributed nuclear matter. In the QMC model it is assumed that each quark feels 
scalar and vector V£ potentials, which are generated by the surrounding nucleons, as 
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well as the confinement potential V(r). Since the typical distance between two nucleons 
around normal nuclear matter density (p = 0.15 fm -3 ) is surely larger than the typical 
size of the nucleon (the radius is about 0.8 fm), the interaction (except for the short-range 
part) between the nucleons should be colour singlet, namely a meson-exchange potential. 
Therefore, this assumption seems appropriate when the baryon density p_B is not high. (For 
high density, see Ref. [[H]].) If we use the mean-field approximation for the meson fields, 



Eq.(|l]) may be rewritten as 

[i 7 • d - (m q - V s q ) - V(r) - ioV v q ]%{r) = 0. (14) 

The potentials generated by the medium are constants because the matter distributes uni- 
formly. As the nucleon is static, the time-derivative operator in the Dirac equation can be 
replaced by the quark energy, — iE. 

By analogy with the procedure applied to the nucleon in QHD [[/]], if we introduce the 
effective quark mass by m* = m q — the Dirac equation Eq.(|T4}) can be rewritten in the 
same form as that in free space, with the mass m* and the energy E — instead of m q 
and E. In other words, the vector interaction has no effect on the nucleon structure except 
for an overall phase in the quark wave function Eq.(f|), which gives a shift in the nucleon 
energy. This fact does not depend on how to choose the confinement potential V(r). Then, 
the nucleon energy (at rest) E^ in the medium is 



E N = M* N (V s q )+Wl (15) 
where the effective nucleon mass depends on only the scalar potential in the medium. 



This important observation was also obtained by non-local field theory model ||18|| , where 
a non-local a-u model containing short distance vertex form factors is used to simulate an 
underlying QCD substructure. 

Although we have discussed the QMC model using the specific model, namely the bag 
model, in our previous works the qualitative features we have found may be correct in 

any model in which the nucleon contains relativistic quarks and the (middle- and long-range) 
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attractive and (short-range) repulsive N-N forces have Lorentz-scalar and vector characters, 
respectively. To confirm this point we will study nuclear matter using the quark models 
discussed in the previous section. 

In Eq.(|T4D the scalar and vector potentials are respectively mediated by the a and u 
mesons. Their mean-field values are introduced by V q = g q a a and V q = g^u, where g q a (g^) 
is the coupling constant of the quark-a (u) meson. (We here ignore the p meson and the 
Coulomb field, for simplicity. For those fields, see Ref. @].) Then, in MFA the effective 
Lagrangian density in the QMC model would be given 

£qmc = #7 • d - M* (<r) - g ulo ojty - ^[(Va) 2 + m 2 a a 2 } + ^[(Vu;) 2 + m> 2 ], (16) 

where ip is the nucleon field, and m a and m u are respectively the masses of the a and u 
mesons — we take m a = 550 MeV and m w = 783 MeV. The u-N coupling constant g w 
is related to the corresponding quark-cj coupling constant as g w = 3g^ [|]. In an uniform 
nuclear matter the derivative terms of the meson fields vanish because the fields do not 
depend on time and position. 

Here we consider the nucleon mass in matter further. The nucleon mass is a function 
of the scalar field. Because the scalar field is small at low density the nucleon mass can be 
expanded in terms of a 

Since the interaction between the quark and the a field is described in terms of the local 
scalar coupling, the derivative of with respect to a is given by 



'dM* N 



SfiJ d? ^ q = -3g^S N (a). (18) 

Here we have defined the quark-scalar charge in the nucleon 5*^(0"), which is itself a function 
of the scalar field, by Eq. (p!8|) . Because Sn is positive, the nucleon mass decreases in matter 
at low density. Note that since the nucleon mass is given by Eq. (|TTD in the HO model the 
derivative turns to 



We furthermore define the scalar-charge ratio, Sn(<t)/ Sn(0), to be C;v(cr) and the er-N 
coupling constant g a : 

C N {a) = S N (a)/S N (0) and g a = 3g*S N (0). (20) 

Using those quantities, we find that the nucleon mass in matter is 

M N = M N - g a a - \g a C' N ($)o 2 + ■■■. (21) 

(In the HO model, we find M N = M N — \g a v — \9oC' N (ti)a 2 + ••-.) In general, Cn is 
a decreasing function because the quark in matter is more relativistic than in free space 
(because of the attractive force due to the a). Thus, C' N (0) takes a negative value. If the 
nucleon were structureless Cn would not depend on the scalar field, that is, Cn would be 
constant (Cn = 1)- Therefore, only the first two terms in the right hand side of Eq . (pT|) 
remain, which is exactly the same as the equation for the effective nucleon mass in the 
original QHD model 0. 



B. Numerical results 

We consider an isosymmetric nuclear matter with Fermi momentum fop, which is given 
by Pb = 2kp/3ir 2 . From the Lagrangian density Eq. fll6|) , the total energy per nucleon E tot 
can be written 

Etot = — ^ / dk^Mfi +k 2 + -^a 2 + ^p B , (22) 
pB^ny J 2p B 2m^ 

where the effective nucleon mass in matter is calculated by the quark model. The value of 
the uj field is determined by baryon number conservation as uo = g LU p B /fn 2 J ||]. 
The scalar mean-field is given by a self-consistency condition 0] 



* 



if* 

dk . N = , (23) 



±9«C N {a) r k r r M N 
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where the derivative of the effective nucleon mass is given by the quark-scalar charge Sn(ct). 
(See Eqs.fllijD and fl20|). For the HO model we should replace CV by §Cjv-) 

In the previous section we have prepared the relativistic constituent quark model and 
the MIT bag model to describe the nucleon. We use m q = 300 MeV in the constituent 
quark model, while the quark is assumed to be massless in the bag model. In the SW model 
we set the radius of the potential to be R = 0.8 fm and detrmine z so as to fit the free 
nucleon mass (= 939 MeV). The parameter (3 is chosen to be (the confining potential 
is purely Lorentz-scalar), 0.25 and 0.5 to examine the effect of Lorentz- vector type confining 
potential. The quark eigenvalue x can be found by solving the boundary condition Eq.([5]) 
at the potential surface. In the HO model, there are two adjustable parameters, c and E s . 
Those are determined so as to fit the free nucleon mass (see Eq. (|TT|) ) and the root-mean- 
square (rms) (charge) radius of the free nucleon: r q {= (r 2 ) 1 / 2 ) = 0.6 fm |2(| (see Appendix). 
The quark energy E is given by the condition Eq. (|10"D . We found that c = 1.591 fm~ 3 , 
E s = 344.7 MeV and E = 641.8 MeV for the free nucleon. In nuclear matter we keep c and 
E s constant and the quark energy E varies depending on the scalar field. 

As in Ref. |3J], in the bag model the bag constant B and the parameter z are fixed to 
reproduce the free nucleon mass under the condition that the nucleon mass be stationary 
under variation of the bag radius R. As in the SW model, we choose the bag radius of 
the free nucleon to be 0.8 fm. (Variations of the quark mass and the radius only lead to 
numerically small changes in the calculated results @.) We found B 1 / 4 = 170.3 MeV for 
massless quarks. In Table | we list up the parameter z, quark eigenvalue and rms radius of 
the free nucleon in the SW and bag models. 

Now we are in a position to determine the coupling constants: g 2 and g 2 are fixed to fit 
the binding energy (—15.7 MeV) at the saturation density (p = 0.15 fm -3 ) for symmetric 
nuclear matter. The coupling constants and some calculated properties for matter are listed 
in Table [TI[ We also present the ratios of the quark eigenvalue, rms radius, axial vector 
coupling constant ((yu), magnetic moment of the proton (/i p ) and quark-scalar charge at 
saturation density to those in free space. (See also Ref. |l{|. For explicit expressions of 
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those quantities, see Appendix.) 

We note that the present quark models (except B(0.25) and B(0.5)) provide the similar 
results in the properties of the in-medium nucleon and nuclear matter. In particular, they 
give values for the nuclear compressibility well within the experimental range (K ~ 200 - 
300 MeV), which is caused by the non-linear dependence of the a field on the nucleon mass. 
In the cases of B(0.25) and B(0.5) the reduction of the nucleon mass in matter is relatively 
small, which is due to the strong Lorentz-vector type confining potential. Furthermore, in 
the HO model the increase of the rms radius at saturation density is about 10%, which is 



large compared with the observed value [21 



We briefly comment on the rms radius, axial vector coupling constant and magnetic 
moment of the nucleon in nuclear matter. In our calculation those quantities are calculated 
using only the (valence) quark wave function and other important contributions from the 
meson (mainly pion) cloud effect etc are not included. In the present quark picture those 
quantities involve integrals over the lower and upper components of the quark wave function. 
Because of the attractive force due to the scalar field the quark becomes more relativistic in 
matter than in free space, that is, the lower component is more enhanced in matter. This 
fact gives the decrease in g^ and the increase in the rms radius and \i v in matter. 

In Fig. [l| we present the total energy per nucleon as a function of nuclear density. Around 
normal nuclear matter density the quark models provide the similar saturation curves, while 
for high density the curves obtained from the models of B(0.25) and B(0.5) are relatively 
lower than those of the others. This corresponds to the fact that a quark model with a 
larger value of (3 gives a lower nuclear compressibility (see Table ||). We show the scalar 
mean-field value and the effective nucleon mass for symmetric nuclear matter in Figs. |2] and 
|3], respectively. From those figures we again see that in the cases of B(0.25) and B(0.5) the 
scalar field is much weaker than those given by the other quark models and the decrease of 
the nucleon mass is fairly small even at high density, which is due to the small scalar field. 

We note that in the SW and bag models with massless quarks the quark-scalar charge 
vanishes in the limit (3 — > 1, which means that the a meson does not couple to the nucleon. 
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(For details, see Appendix.) This fact implies that as (3 is larger the cr-N coupling is weaker 
in matter. Thus, we may conclude that qualitatively a large mixture of the Lorentz-vector 
type confining potential leads to a weak scalar mean-field and a large effective nucleon mass 
in nuclear matter. It is well known that within MFA a smaller effective nucleon mass (and 
hence a stronger scalar field) in matter is favorable to fit various properties of finite nuclei [|7| . 
Therefore, the confining potential including a strong Lorentz-vector type one may not be 
suitable for describing a nuclear system. 



IV. TRANSFORMATION FROM QMC TO QHD 

The main difference between the QMC and QHD models lies in the dependence of the 
scalar field on the nucleon mass in matter. By performing a re-definition of the scalar field, 
the QMC Lagrangian density at hadronic level can be cast into a form silimar to a QHD-type 
mean-field model, in which the nucleon mass depends on the scalar field linearly, with self- 



interactions of the scalar field [[^,|I(J. Here we will briefly review the field re-definition and 
discuss, as an example, the derivative (scalar) coupling model of QHD (i.e., the Zimanyi- 
Moszkowski (ZM) model) ||23|| . Then, we will calculate non-linear scalar potentials generated 
from the QMC model and compare with the potentials phenomenologically determined by 
fitting the properties of finite nuclei and nuclear matter. 

A. Field re-definition 

In the QMC model the nucleon mass in matter is given by a function of the scalar field 

a 

M* N , QMC = M* N>QMC (a), (24) 
while in QHD the mass depends on a scalar field <fr linearly, 

M^ QHD = M N - g 0< j>. (25) 
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This fact suggests a re-definition of the scalar field in QMC: 

g <P(a) = M N -M* Nmw (a) } (26) 
where go is a constant chosen so as to normalize the new scalar field <fi in the limit 0^0: 

0(a) = a + 0(a 2 ). (27) 

Thus, g is given by 

Comparing with Eqs.flTBD and (|20"D, we find go = g a for the SW and bag models and go = \g a 



for the HO model. 

The contribution of the scalar field to the total energy, E sc i, is now rewritten in terms of 
the new field <fi. From Eq.(|16|) we find 

E scl = \jdr [(W) 2 + mla 2 ] = J df [h{<P)\V<pf + U s (<f>)\, (29) 

where U s describes the self-interactions of the scalar field 

U s (cj>) = l -mla{<i>) 2 , (30) 

and 

Note that in uniformly distributed nuclear matter the first derivative term in E sc i does not 
contribute to the total energy. (The effect of this term on the properties of finite nuclei was 
studied in Ref. [Kj.) 



The QMC model is now re-formulated in terms of the new scalr field 0, and it is of the 
same form as a QHD model with the non-linear scalar potential U s (4>), which is generated 
from the effect of the internal structure of the nucleon, and with the coupling h{4>) to the 
gradient of the scalar field. However, we should keep in mind that this re-formulation is 
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performed only at hadronic level and that the non-linear potential U s is simply determined 
so as to simulate the effective nucleon mass given by the quark model in matter. Thus, only 
on the energy of nuclear matter the QHD-type mean-field model derived here is identical to 
the QMC model. We will again discuss this point in the last section. 
The standard form of the non-linear scalar potential is ||[| 

U s {<P) = \ml^ + ^ + ^\ (32) 



In more sophisticated version of QHD p2| , inspired by modern methods of effective field 



theory, many other terms of meson-meson and meson-nucleon couplings are considered. 
However, rather than attempting to compete with those models, in this paper we concentrate 
only on the parameters k and A to make our discussion simple. It is well known that the non- 
linear scalar potential Eq. (|32]) is practically indispensable to reproduce the bulk properties 
of finite nuclei and nuclear matter ||j9l,p2]-p5|| . 

If the effective nucleon mass could be determined experimentally as a function of the 
scalar field, this procedure leading to Eq. (|3~0p would allow to predict a part of the non-linear 
potential, which is originated by the quark substructure of the nucleon in nuclear matter. 
Inversely, if the non-linear potential obtained by fitting various experimental data could be 
separated into a portion due to the effect of the internal structure of the nucleon and the 
inherent self-interaction part of the scalar field, we could get significant information on the 
quark structure of the nucleon in matter by using the inverse field transformation, 



= JL_=^Z. (33) 

It is potentially of very interesting. However, it may be quite difficult. 

As an example, we here consider the derivative coupling model of QHD, namely the 
Zimanyi-Moszkowski (ZM) model ||23|| . By re-definig the scalar field the ZM model can be 
exactly transformed to a QHD-type model with a non-linear potential. 

In the ZM model the interaction term between the a and the nucleon is given by 



£fnf = ■ fy, (34) 
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in addition to a similar coupling to the vector field. Then, after rescaling the nucleon field, 
the Lagrangian density turns to be the same form as that of the original QHD, with the 
effective nucleon mass 

instead of the mass depending on the scalar field linearly. Thus, the model involves higher 
order couplings between the a and the nucleon, which provides a good nuclear compressibility 
around 230 MeV and M* N /M N ~ 0.85 g|. 
Let us define a new scalar field by 

a) = M N - M* (cx), (36) 



(37) 



where the mass Mj^ zmW) ls given by Eq . (|35"D . Then, we easily find 



l + g'a 

with g' = g a /Mjsf. Inversely, we also find 



= j^. (38) 



The non-linear potential is thus given by Eq. (|30| 

2 2 V 1- ^ 

which can be expanded up to order 0(4> A ) 



2 



= »(0) 2 = X ( T^72) > (39) 



17.(0 * l -ml<t? + g„ 3 + \gl { (40) 



2 aT " \M N ) T T a \M N . 
(The shape of the potential Eq.([3U|) is shown in Fig.|5].) Hence, the ZM model is identical to 
the QHD-type model with the non- linear scalar potential expressed by Eq . (|39|) . 

B. Non-linear scalar potentials generated from QMC 

In general the nucleon mass in matter may be given by a complicated function of the 
scalar field. However, in QMC the mass can be parametrized by a simple expression up to 
0{gl) Hi (see Fig. § 
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^L^l-ay + by 2 , (41) 

with a dimensionless scale y = g a a/M^ and two (dimensionless) parameters a and b. From 
the discussion around Eq.fl2lD we can expect that a ~ 1 for the SW and bag models while 
a ~ 2/3 for the HO model. This parametrization is accurate up to ~ 4p - Using Eq.(|41"|) the 
calculated nucleon mass in nuclear matter can be fitted with the method of least squares. 
We show the parameters a and b for the various quark models in Table [TTT| . We find that 
the present quark models (except B(0.25) and B(0.5)) give b ~ 0.2 — 0.5. 

Now we shall re-define the scalar field using Eq. (|26|) . We find go = ag a and 



(f)(a) = a- da 2 , (42) 
with d = bg a jaM^. Eq . (|42|) can be solved exactly, and we obtain 



1 - VI - 4d0 

<4>) = ^ ' (43) 

which satisfies the condition a —>■ in the limit 0^0. The non-linear potential is then 
calculated 

U M = f , (44) 

and we can expand C/ s as 

W) = !^ + 9 „ r (g / + 5 9 , r2 + 0{g% (45) 

where r = b/a. 

Comparing with Eq.(^) we can find the parameters k and A in the non- linear scalar 



potential generated by QMC They are listed in Table |TV|. The quark models we have 
adopted here lead to the non-linear potential with k ~ 20 — 40 (fm _1 ) and A ~ 80 — 400, 
except for the models of B(0.25) and B(0.5). (The models B(0.25) and B(0.5) give quite 
large values of k and A, in addition to the large effective nucleon mass and the small scalar 
field in matter, which seems impractical.) From the discussion in the subsection A| we 
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see that both a and b are positive in QMC because a and b are respectively proportional 
to SV(0) and —C' N (0). Hence, from Eq.fl4"5|) we can expect that the effect of the quark 
substructure of the nucleon in matter is cast into a non-linear potential with positive k and 
positive A in the QHD-type mean-field model. 

In Table [V], for comparison we show the parameters k and A, which are phenomenologi- 
cally determined in various relativistic mean-field (RMF) models. In such models k and A 
are sometimes taken to be negative. Furthermore, in Figs. |] and ||] we respectively illustrate 
the shapes of the non-linear scalar potentials calculated from the QMC model and those 
in the RMF models. (In Fig. |] we did not show the results of B(0.25) and B(0.5).) From 
Fig. f| we can see that the various quark models lead to the similar non-linear potentials, 
despite of the big difference in the confinement mechanism. It may imply that we cannot 
discriminate the quark models for the nucleon only from the point of view of the energy of 
nuclear matter. On the other hand, the non-linear potentials in the RMF models show quite 
different behaviors in Fig. |5|. In the region of <p > the potentials in the G2, NLB and ZM 
models are close to those produced by QMC. 

We here comment on the coupling h(4>) to the gradient of the scalar field. The QMC 
model gives 

ftW = 1+s -'te) +to (^) 2 + °(^ < 46) 

where g s \ = 2g a r and g S 2 = 6g%r 2 . We find that the present model suggests g s i ~ 4 — 8 
and g S 2 ~ 20 — 100, while the ZM model gives g s \ ~ 30 and g S 2 ~ 580. In the RMF model, 
for example, the Gl and G2 models lead to g s % ~ 9, which is consistent with the values 



obtained in QMC. This coupling h(4>) affects the surface properties of finite nuclei [JT0 



V. DISCUSSION AND CONCLUSION 

We have first studied the MIT bag model and the relativistic constituent quark model 
for the nucleon, where a square well and harmonic oscillator potentials are used to confine 
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the quarks. We have considered not only the Lorentz-scalar type confining potential but 
also the Lorentz- vector type one in order to study the role of the vector potential in nuclear 
matter. Then, we have calculated the properties of the nucleon and nuclear matter using 
QMC with those quark models. As a mixture of the Lorentz- vector type confining potential 
is larger the effective nucleon mass in matter is larger and the strength of the scalar field 
is weaker. In particular, it is preferable in the bag model that the Lorentz structure of the 
confinement is purely scalar (or including a very weak Lorentz- vector type potential). 

Next, we have performed a re-definition of the scalar field in matter and transformed the 
QMC model to a QHD-type model with a non-linear scalar potential. We then compared our 
results with the potentials which are determined so as to fit the properties of finite nuclei and 
nuclear matter in the RMF models. The QMC model provides the parameters k ~ 20 — 40 
(fm _1 ) and A ~ 80 — 400 for the non-linear scalar potential. The shapes of the potentials 
generated from QMC are thus very close to one another, despite of the different confinement 
mechanism in the quark models. On the other hand, the parameters phenomenologically 
determined in the RMF models take various values and the shapes of the potentials are quite 
different. In general, the phenomeno logical potential may consist of a part, which is caused 
by the quark substructure of the nucleon, and the inherent self- couplings of the scalar field. 
Therefore, if the part due to the internal structure of the nucleon could be known by nuclear 
experiments, we could get significant information on the quark structure of the in-medium 
nucleon. 

We should emphasize that the field re-definition performed here enables us to make a 
QHD-type effective Lagrangian from the QMC model. However, the created Lagrangian is 
identical to the QMC model only on the energy of a nuclear system. In QMC the quark 
wave function in the nucleon is modified by the nuclear environment self-consistently at 
each nuclear density. We can know it explicitly. It is the big difference between the QHD- 
type and QMC models, and is of great advantage to the QMC model. Using the change of 
the quark wave function in the medium provided by QMC we can explain many intriguing 
nuclear phenomena ||, for example, the nuclear EMC effect |26], charge symmetry breaking 
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in nuclei (the Okamoto-Nolen-Schiffer anomaly) p7j , change of electromagnetic form factors 



of the bound nucleon |55[], longitudinal response functions in electron scattering |29[] etc. 

We have one comment on the modified quark-meson coupling (MQMC) model proposed 
by Jin and Jennings The original QMC model |3|-|5| generates a (relatively) small 

scalar field in nuclear matter. On the other hand, it is well known that the MQMC model 
can generate the scalar field, which is close to that in QHD (see also Ref. |10|). However, 
unfortunately the increase of the nucleon size at po is about 40% of the free-space value, 



which is apparently too large |2l||28]]. Although it must be true that the bag constant 
decreases as the density grows up, the reduction of the bag constant around normal nuclear 
matter density may not be large. For example, the Dyson- Schwinger approach at finite 
chemical potential |31| suggested that the bag constant is reduced very little at p , while it 
decreases very rapidly near the phase transition point. 

It is possible to extend the present QMC model to a model which includes quark degrees 
of freedom in the mesons, as well as in the nucleons — we call it QMC-II P,|32"f . The QMC- 
II model can provide a lot of effective coupling terms among the meson fields because the 
mesons have structure. In particular, the QMC-II Lagrangian automatically offers non-linear 
terms with respect to the meson fields. Those interaction terms may correspond to the terms 
of meson-meson and meson-nucleon couplings considered in the new version of QHD |22 



or 



higher order terms appearing in the chiral effective lagrangian for nuclear matter [Q. It 
will be very interesting to construct a QMC model, in which both the quark substructure 
of hadrons and non-linear potentials due to inherent self-couplings of the meson fields are 
involved (see, for example, Ref. [[12]]), and explore the connection between various coupling 
strengths found empirically in the phenomenological models and those predicted in the QMC 
model. 
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APPENDIX 



In the relativistic quark models we adopt here the lowest energy state for the confined 
quark is given by Eq. (Q) . In the SW model the upper and lower components of the quark wave 
function inside the well are respectively given by the spherical bessel functions: f(r) = jo{kr) 



and g(r) = J (a — A)/(a + X)ji(kr), where a = RE = y/x 2 + A 2 , A = m q R and k = x/R. 
The normalization constant is then given by 

N 2 = (47) 

where D sw = 2a (a - ^l- (3 2 + (3\) + (3 2 . 

Using those quantities the quark-scalar charge Sn is calculated 



SN(m q ) = dr ip q i[) q , 
Jv N 



2A(a - ^JT^W + /3\) + a^l^ 1 



(48) 



Dsw 

We note that vanishes if m q = MeV and (5 — > 1, which may be accidental in the SW 
and bag models. In this case there is no coupling between the a meson and the nucleon, 
which implies that the nuclear matter cannot be described by this model. 
The root-mean-square (rms) (charge) radius of the nucleon is given by 

r\ = (r 2 ) = / dr r^ q = R 2 -%^, (49) 
>v N 2>x 2 D S w 



where 



N rms = a[2x 2 (a -\Jl- /3 2 ) + 2a(2 + /3) - 3] - 3A 
+ 2aA 



A(l + 2/5)-a: 2 v /l-/5 2 - 2 



j 3 



x 2 --)P-2 



(50) 



The axial vector coupling constant (gX) and the magnetic moment (/x) of the nucleon are 
also calculated analytically. In the present quark model qa is given by g& = (5/3) x g° A) 
where the factor (5/3) is the expectation value of the spin-isospin operator with respect to 
the nucleon state. Then, g° A is 
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_ 2a 2 + AaX - 3X^/1 - f3 2 + 2(3X(a + 2X) 
9a — TFj • (51) 

The magnetic moment is given by ft = Y.i^i^iQi^ where Qi is the quark charge operator 
and the sum i runs over all quarks in the nucleon. Then, /ij is 

R Aa + 2X- 3y/l - (3 2 + 2(3{a + 2A) 
^ = fi" n • ( 52) 

O L>sw 

In the MIT bag model, those quantities are given by the same expressions as in the SW 
model. 

In the HO model, the quark wave function is again given by Eq.(^j). The upper and 
lower components are respectively given by gaussian functions: f(r) = e~ r ' 2 ^ 2r o and g(r) = 
(r/£r 2 )e _T ' 2/ ' 2r o, where r$ = (c^)^ 1 / 4 with £ = E + m q . Then, the normalization constant is 

^ 2 = 4 • (53) 

y>H L> HO 

where D HO = 2(£r ) 2 + 3. 

Using the quark wave function, the quark-scalar charge is calculated 

2(£r ) 2 - 3 

SN[m q ) = = • (54) 

L>HO 

We note that Sn approaches 1/3 when m q = MeV (and (3 = 1). Similarly, g A and \ii are 
respectively given by 

A = ^f^, (55, 
/* = jf^- (56) 

DUO 

When the nucleon is embedded in a nuclear medium, those quantities can be expressed 
by the same forms with the effective quark mass m* (see below Eq.(^)), instead of the free 
quark mass m q . 
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TABLES 



TABLE I. 


Parameter 


z, quark eigenvalue 


x and rms 


radius (r q ). The ba 


g model is 


denoted by 


B. The value of /3 is shown in the parentheses. 










SW(0) 


SW(0.25) 


SW(0.5) 


B(0) 


B(0.25) 


B(0.5) 


z 


4.396 


4.792 


5.164 


3.273 


3.974 


4.637 


X 


2.449 


2.596 


2.732 


2.043 


2.276 


2.498 


r q (fm) 


0.5415 


0.5264 


0.5137 


0.5832 


0.5721 


0.5593 



TABLE II. Coupling constants and calculated properties of the nucleon and symmetric nuclear 
matter at normal nuclear density. (The quantities with * are those quantities calculated at p$.) 
The nuclear compressibility K is quoted in MeV. The last row is for QHD J?J. 





9l 


9l 


M^/M N 


K 


x*/x 


r\jr q 


9a/ 9A 




Sjv/Sn 


SW(0) 


84.41 


104.1 


0.725 


329 


0.937 


1.03 


0.936 


1.11 


0.886 


SW(0.25) 


74.75 


83.16 


0.769 


307 


0.956 


1.03 


0.919 


1.11 


0.850 


SW(0.5) 


66.61 


65.21 


0.807 


287 


0.971 


1.03 


0.903 


1.10 


0.811 


HO(l.O) 


146.5 


64.51 


0.805 


309 




1.10 


0.891 


1.23 


0.812 


B(0) 


67.55 


66.10 


0.805 


278 


0.837 


1.02 


0.909 


1.08 


0.814 


B(0.25) 


54.65 


35.49 


0.869 


235 


0.900 


1.03 


0.877 


1.06 


0.711 


B(0.5) 


47.99 


10.92 


0.921 


169 


0.941 


1.02 


0.841 


1.05 


0.541 


QHD 


91.64 


136.2 


0.556 


540 
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TABLE III. Parameters a and b. Note that a ~ 1 for the SW and bag models, while a ~ 2/3 
for the HO model. 





SW(0) 


SW(0.25) 


SW(0.5) 


HO(1.0) 


B(0) 


B(0.25) 


B(0.5) 


a 


1.008 


1.006 


1.017 


0.6868 


0.9982 


1.005 


0.9933 


b 


0.2146 


0.3229 


0.4965 


0.2452 


0.4346 


0.9666 


2.206 



TABLE IV. Parameters n and A generated by the QMC model. 





SW(0) 


SW(0.25) 


SW(0.5) 


HO(1.0) 


B(0) 


B(0.25) 


B(0.5) 


n (fm 1 ) 


19.16 


27.20 


39.03 


42.31 


35.05 


69.65 


150.7 


A 


78.78 


158.7 


326.8 


384.1 


263.5 


1041 


4872 



TABLE V. Parameters k and A for the RMF models. 



ZM j23| 



NL1 



TM1 |5[ NLB p 



Gl 1221 



G2 1231 



k (fm- 1 ; 

A 



74.8 
720 



-24.3 
-218 



-14.5 
3.71 



4.06 
10.0 



29.4 
-279 



46.0 
19.7 
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FIGURES 

120 i , , , 




Pb/Po 

FIG. 1. Energy per nucleon for symmetric nuclear matter. The saturation curves for the SW 
model and the bag model are respectively shown by solid and dotted curves. The dot-dashed curve 
presents the result of the HO model. The curves with squares (triangles) are for the results with 
(3 = 0.25(0.5), while the curves without any marks are for (3 = 0. 
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Pb/Po 

FIG. 2. Scalar mean-field value. The curves are labeled as in Fig. |l[ 
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FIG. 4. Non-linear scalar potentials generated from the QMC model. The solid curve shows 

2 

U s = ^f-cf 2 . The dashed curve with open (solid) circles is for the SW model with f3 = 0.25(0.5), 
while the dashed one without any marks is for the SW model with (3 = 0. The result of the B(0) 
model is shown by the dotted curve. The dot-dashed curve is for the HO model. 
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2 

FIG. 5. Non-linear scalar potentials in the RMF models. The solid curve shows U s = ^-cf 2 . 



The dashed curve with open (solid) circles is for the G2 (NLB) model 22], while the dashed 



one without any marks is for the Gl model [22|. The dashed curve with crosses is for the ZM 



model [23]. The potentials in the TM1 [25] and NL1 24] models are respectively shown by the 
dotted and dot-dashed curves. 
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